Mottness on a triangular lattice 
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We study the physics on the paramagnetic side of the phase diagram of the cobaltates, N a x CoOi, 
with an implementation of cellular dynamical mean field theory (CDMFT) with the non-crossing 
approximation (NCA) for the one-band Hubbard model on a triangular lattice. At low doping we 
find that the low energy physics is dominated by a quasi-dispersionless band. At half-filling, we find 
a metal-insulator transition at U c = 5.6 =b 0.152 which depends weakly on the cluster size. The onset 
of the metallic state occurs through the growth of a coherence peak at the chemical potential. Away 
from half filling, in the electron-doped regime, the system is metallic with a large, continuous Fermi 
surface as seen experimentally. Upon hole doping, a quasi non-dispersing band emerges at the top 
of the lower Hubbard band and controls the low-energy physics. This band is a clear signature of 
non-Fermi liquid behavior and cannot be captured by any weakly coupled approach. This quasi 
non-dispersive band, which persists in a certain range of dopings, has been observed experimentally. 
We also investigate the pseudogap phenomenon in the context of a triangular lattice and we propose 
a new framework for discussing the pseudogap phenomena in general. This framework involves 
a momentum-dependent characterization of the low-energy physics and links the appearance of 
the pseudogap to a reconstruction of the Fermi surface without invoking any long range order or 
symmetry breaking. Within this framework we predict the existence of a pseudogap for the two 
dimensional Hubbard model on a triangular lattice in the weakly hole-doped regime. 



I. INTRODUCTION 

Charge carriers in the cobaltates, Na x Co02, are lo- 
cated in two dimensional C0O2 layers separated by insu- 
lating layers of iVa + ions which act as electron donors. Their 
structure is a triangular net of edge-sharing oxygen octahedra 
with the Co atoms occupying the center and the Na atoms 
playing the role of electron donors. The octahedral symmetry 
around the Co ions results in a splitting of the d-orbitals in 
two e g and three lower lying ti g orbitals. The trigonal dis- 
tortion of the C0O2 layers further splits the t<i Q orbitals 
into one a\ g and two lower e' g Singh 1 . The valence of the 
cobalt ate ions is Co A ~ x which means that the Fermi sur- 
face will lie in the a\ g orbital which will range from half 
to fully filled. Consequently, the cobaltates constitute a 
realization of strongly correlated electron physics on a tri- 
angular lattice. Across their phase diagram they exhibit 
a wide range of behavior 2 ranging from a paramagnetic 
Fermi liquid at low Na concentration x, a strange Fermi 
liquid with Curie Weiss magnetic susceptibility for high 
x and a singular insulating state at x = 0.5. While the 
paramagnetic metal exhibits some properties akin to that 
of a Fermi liquid, the cobaltates still remain strongly cor- 
related systems. For example, experiment- and theory^ 
place the hopping matrix element and the on-site repul- 
sion at t = 0.2eV and U = 4eV, respectively. 

Various ARPES studies have been performec&&£, 
which suggest a Fermi surface which consists only of a 
large a\ g hole pocket with the e' g orbitals lying under the 
Fermi surface. This is in contrast with LDA calculations 
Singh— which suggest the existence of peripheral e' g hole 
pockets. In an effort to resolve this discrepancy, several 
multi-band DMFT studies have been performed^. Fi- 



nally a few CDMFT calculations Kyung 10 , Lee et al.— 
address mostly the Mott transition on a triangular lattice 
or compare different impurity solvers. 

Motivated by the cobaltates, we address in this paper, 
the properties of strong electron correlation on a trian- 
gular lattice. Of particular interest will be the nature of 
the Mott transition at half-filling on such a lattice. A 
triangular lattice offers an ideal playground for exploring 
the Mott transition as a result of the inherent magnetic 
frustration that is present. We find that a critical value of 
U = b.7t separates paramagnetic insulating and metallic 
phase. Away from half-filling we find a metallic phase 
with a large Fermi surface as is observed experimentally. 

This paper is organized in three main sections. In sec- 
tion [Til we give an overview of the computational scheme, 
the cluster dynamical mean field theory with the non- 
crossing approximation as the impurity solver in the con- 
text of the one-band Hubbard model. In section [Till we 
discuss the the two main issues related to the consistency 
of the method: the proper periodization procedure to ob- 
tain physically correct lattice quantities and the cluster 
size dependence of the the results. We show that, for a 
given small cluster size, the method breaks down at cer- 
tain filling values and we argue that the cluster size in- 
dependence should be the ultimate consistency criterion. 
In section [IV] we present the results of the simulation. 
In IIV At we show the existence of a quasi-dispersionless 
low energy band, which is signature of strong correlations 
and incompatible with Fermi-liquid physics. In IIV B[ we 
discuss the Mott transition on the triangular lattice. Fi- 
nally, in IIV C\ we argue for the existence of a pseudogap 
at low hole doping. 
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II. DESCRIPTION OF THE METHOD 

We start with the one-band Hubbard model, 

h = -t c L c j* + c - c - + u w 

<r,(hO) i 

where t is the matrix element for hoppings between near- 
est neighbor sites, and U is the on-site repulsive 
interaction. We assume that the single-band Hubbard 
model on a triangular lattice captures the main features 
of strongly correlated physics in the presence of magnetic 
frustration. We restrict our study to the paramagnetic 
state 

(nil) = <"iT> = »/2 (2) 

which is consistent with the experimental observations in 
cobaltates for x < 0.5. 

As a computational tool, we use in our investigation a 
real-space cluster generalization of dynamical mean field 
theory (DMFT).— The DMFT has been a very successful 
tool in investigating many aspects of strongly correlated 
systems. In this method one single site is treated as an 
impurity embedded in an effective bath consisting of the 
rest of the sites the properties of which are captured by 
the hybridization function. It is exact in infinite dimen- 
sions or more precisely in infinite coordination number z 
and it can successfully describe the antiferromagnetic or- 
der. However for many applications it is necessary for the 
short range (few lattice site) correlation to be described 
accurately. In the cluster DMFT method (CDMFT)^ a 
cluster extending in a small number of sites is treated as 
the impurity and therefore the local (cluster) degrees of 
freedom are treated exactly. The rest of the lattice, the 
bath, is described by a multi-component hybridization 
function. 

All cluster-DMFT-based algorithms contain the follow- 
ing major components. 

• An impurity solver, which evaluates the cluster 
Greens function from the hybridization function. 

• A self consistency condition which expresses the 
hybridization function with respect to the cluster 
Greens function. 

• A periodization procedure which connects the lat- 
tice quantities with the cluster quantities. 

A. Impurity solver 

The impurity solver evaluates the cluster Green's func- 
tion, given the "external" hybridization function. Vari- 
ous impurity solvers have been proposed in the literature 
such as exact diagonalization and quantum Monte- Carlo. 
However those can only be implemented in imaginary 
time and an analytic continuation is required to obtain 
real time properties. A real time impurity solver is the 



non-crossing approximation (NCA), which is a first or- 
der perturbation theory with respect to the hybridization 
function. It has the advantage of being very fast and rel- 
atively easy to implement. The NCA has been a valuable 
tool for extracting the physics of the Anderson impurity 
models. The NCA equations can be obtained by using 
the slave boson method 14 and they can be expressed with 
respect to the pseudo-particle Resolvents G mn and their 
self energies £ mn where m, n are indices representing the 
eigenstates of the cluster. The NCA equations, which 
are used to evaluate the updated resolvent self-energies 
along with the cluster Green function G^ v for a given 
hybridization function A uu , are 1 ^: 

m'n' \iv 
m'n 1 uv 

J de/(-0A M „(0G m '„' (iu - (3) 

g„„m = -i £ F;'"(F-' m )*x 

mnm'n' 

-G m . n ,(Z-iu)G nm (Z)] (4) 
Q = - /^e-«£G rorn (0 

m 

The Greek indices correspond to cluster degrees of free- 
dom (site and spin), F™ m = (m\c u \m f ) are the matrix 
elements of the destruction operator and /(£) is the Fermi 
function. The resolvents can be obtained from the self 
energy using the Dyson equation, 

G mn {iw) = (iu-E-X-Yi)^ n (5) 

where E is the diagonal matrix of the clusters eigen- 
energies and A, an artifact of the slave-boson approach, 
is chosen for convenience. 

In the paramagnetic, disordered state, both the spin 
and the irreducible representation of the geometrical 
symmetry group are good quantum numbers and can be 
used to label the cluster eigenstates, 

\m) = \N,S\S z ,r,r m ,E) (6) 

where r and r m denote the irreducible representation and 
it's row, respectively and the rest follow standard nota- 
tion. All but the energy E are good quantum numbers 
and they cannot be affected by the bath. Furthermore 
all resolvents with the same N, 5 2 , r but different S z and 
r m are the equal to one another. 
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B. Self consistency condition 

The self consistency condition generates a new hy- 
bridization function, i.e., a new effective environment 
for the impurity cluster, starting from a given cluster 
Green function G^ Vl and taking into account the geome- 
try of the lattice and the non-interacting hoping matrix. 
Within the CDMFT scheme, the self consistency condi- 
tion reads 



^(M- 1 - E(K))- 1 = (M _1 - A - To)" 



(7) 



K 



where the matrix E(K) is the Fourier transform of the 
inter-cluster hoping matrix X/j, with / and J being in- 
dices that label the clusters inside the superlattice, and 
K is the superlattice momentum. The cluster cumulant 
M a p can be expressed in terms of the intra-cluster hoping 
matrix Tq = Tjj as 



M~ 



G~ 



A + T . 



(8) 



C. Lattice periodization 



Once the cluster quantities such as the cumulant or the 
self energy have been obtained, the corresponding lattice 
quantities need to be reconstructed D. et al.— . A good 
estimate of these lattice quantities can be obtained by 
averaging over all the possible ways in which a lattice 
can be covered with clusters of a given type. For each 
type of cluster and each lattice there are a finite number, 
Nsj of different realizations of the superlattice, which 
are related to one another by a symmetry operation (ro- 
tation, translation or both). Explicitly, a lattice quantity 
Xiatt{x-a — x /3) 5 which may be either the cumulant of the 
self energy, can be extracted from the corresponding set 
of cluster components X(x a ,x^) as 



Xlatt^o 



(*!%]) (9) 



where X SL represent the quantity X for a certain ref- 
erence superlattice, S is a symmetry operation relating 
different equivalent superlattices to the reference super- 
lattice and S[x a ] is the new position of site x a after ap- 
plying the symmetry operation. If the positions 5[x a ] 
and S[xp] do not belong to the same cell of the refer- 
ence superlattice X SL (S[x a ], S[xp]) vanishes, otherwise 
it is given by the corresponding cluster component. The 
momentum dependence is obtained by a simple Fourier 
transform, 

X P (k) = J2 X iatt(** ~ x^)e^ x «- x ^- k (10) 

P 

where the index P signifies that the quantity X(k) was 
obtained by applying the periodization procedure. 



In this study, we focus on the cumulant periodization 
scheme X = M, in which the lattice Green function is 
given by 



M-\u;,k)-e(k) 



(11) 



We also discuss briefly the implications of using the self 
energy periodization scheme. Note that the cluster self 
energy is related to the cluster cumulant through the 
matrix equation 



E = (lj + /i) 1 - AT 



(12) 



and a similar scalar equation holds for the corresponding 
lattice quantities, £(k) = {uj + /i) — M(k) _1 . In this 
scheme the lattice Green function is given by 



G(cj,k) 



uj + \± — e(k) — Sp(cj.k) 
1 

(M- 1 ) p (^,k)-e(k)' 



(13) 



where the index P implies that M _1 is periodized. 

In this study we are going to use two types of clus- 
ters: a triangular 3-site cluster and a rhombic 4-site clus- 
ter. Because of paramagnetism and also the geometrical 
symmetry, the triangular cluster has 31 independent re- 
solvents and 2 independent cluster quantities, whereas 
the rhombic cluster has 309 and 5 respectively. Single 
site and two-site clusters have also been considered but 
convergence is possible only at higher temperatures. The 
respective superlattices are shown in Fig. O 

For the triangular cluster, there are only two indepen- 
dent components, a local Xo and the nearest neighbor 
one X\. The corresponding periodization is 



X tri (k) =X + 2Xia(k) 



where a(k) = ^Xli=i cos ^ an d &i = k x , k<i — 7} > 
^ky, ks = — \k x — ^-ky. For the rhombic cluster there 
are 5 independent components: 2 local ones Xo and Xq 
corresponding to the site with 3 and 2 neighbors inside 
the cluster respectively, two nearest neighbors X\ and 
X{ corresponding to one of the sides of the cluster and 
the diagonal link respectively and one next-to nearest 
neighbor component X2 along the long diagonal of the 
cluster. The periodization is: 



X r ho(k) 



Xx 



X' 2 



2X 1 + -X[ a(k) 



X 2 



6(k) 



where b(k) = | Y^l=i cos — fe+i) with k± = k\. Since 
a uniform paramagnetic phase is assumed, the choice of 
the cluster is expected to have a relatively small impact 
on the physical quantities for regimes characterized by 
short correlation lengths. 
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Figure 1: Types of clusters used in the simulation: (a) triangle 
(b) rhomboid. 




Figure 2: Triangular (left) and Rhomboid (right) cluster's 
superlattice. The black circles represent the superlattice's 
sites and the dashed arrows its unit vectors. 




— I — 
new A 
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Figure 3: Flow diagram for the CDMFT+NCA method. A 
is the bath function, E is the resolvent self energy and G the 
cluster Green function and the cumulant is M = G -1 +A+Tb. 
The flow starts at the Guesses of A (usually with a Gaussian 
imaginary part) and E (negative imaginary constant). The 
impurity solver evaluates a new E from A and the old E and 
also the cluster Green function. The self consistency condition 
takes as its input the cluster Green function and the old A 
to return an updated A. The impurity solver and self consis- 
tency condition iterate until convergence is reached. Then the 
cumulant M is evaluated and periodized. Physical quantities 
can be obtained from M (for example, the spectral function) 
and E. 



III. NUMERICAL SCHEME: CONSISTENCY, 
OPTIMIZATION AND LIMITATIONS 

The key feature that makes a DMFT-type treatment 
applicable is the locality of the correlated physics. In infi- 
nite dimensions, the correlations are purely local and can 
be described by a momentum-independent self-energy. In 
finite dimensions, the basic assumption is that the corre- 
lations are short ranged and can be captured by a cluster 
extension of DMFT. The size of the cluster that would 
properly capture the physics is determined by the range 
of the relevant correlations and cannot be known a pri- 
ori. Therefore, consistency checks are a necessary com- 
ponent of any cluster DMFT treatment. In this section, 
we show that for the two-dimensional Hubbard model on 
a triangular lattice: A) The self-energy is not a short- 
range quantity in the Mott insulating phase and nearby 
and therefore should not be extracted from the cluster 
components. Instead, the renormalized two-point cumu- 
lant satisfies the locality requirement and can be used 
for re-constructing the lattice quantities. B) A cluster 
scheme does not work equally well for all doping val- 
ues. In particular, for certain doping levels commensu- 
rate with the cluster size the scheme predicts spurious 
"insulating" states. We argue that a comparison between 
results obtained using clusters of different sizes is cru- 
cial. The ultimate consistency criterion is the invariance 
of results to an increase of the cluster. 

To solve for the cluster quantities and the resolvents 
self consistently, we start from an initial guess for the 
imaginary part of the resolvent self energies and the hy- 
bridization function. The real part was obtained through 
the Krammers-Kronig relationships. One possibility is to 
start at high temperatures (T w 0.3t) where the method 
converges very easily (a constant ^E mn and a Gaussian 
^A^ u is enough), and then "cool down" progressively us- 
ing in every step the solution of the previous step. Usu- 
ally a two step process suffices. Once the initial guess 
is obtained, the NCA equations, (0 and ((4]), along with 
Eq. ([5]) are used to update E mn and evaluate the cluster 
Green function G^ v . From the latter and the self con- 
sistency condition, Eq. 0, the hybridization function 
A^v is updated. The process iterates until convergence is 
reached. The lattice Green function and the correspond- 
ing spectral function are obtained by the periodized cu- 
mulant. A flow diagram of the process is shown in Figure 
03 From the spectral function, a variety of two-particle 
properties can be obtained. 



A. Cumulant versus self-energy periodization 

Our first task is to test the accuracy of the periodiza- 
tion procedure and identify the quantity most suitable 
to be used in the periodization scheme. Note that the 
implicit physical assumption behind periodizing a cer- 
tain quantity X(x a , x^) is the short-range nature of that 
quantity. A long-range quantity cannot be properly ap- 
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Triangle, n=0.99 



Rhomb, n=0.99 




n=0.990 



n=1.010 




n=1.160 




■ Cumulant 
- Self Energy 



n=1 .492 




Figure 4: Comparison of the local cluster spectral function, 
An, with the one obtained through cumulant Mp and self- 
energy Hp periodization, for the triangular (left) and rhombic 
cluster (right) at U = 12t and T = OAt. The local cluster 
spectral function agrees with that from Mp. The one obtained 
from Hp predicts unphysical states in the middle of the Mott 
gap. 
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n=1.145 




n=1.008 
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Figure 5: Comparison of the density of states for the cumulant 
and self energy periodization for the triangular cluster and 
U=12t, T=0.1t. 



Figure 6: Comparison of the density of states for the cumulant 
and self energy periodization for rhombic cluster at U=12t, 
T=0.1t. 
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Figure 7: The density of states for both clusters for (a) n 
1.22 and (b) n = 1.33 



odization procedure. This signals that the cumulant is a 
short-range quantity in the vicinity of the Mott transi- 
tion, while the self-energy is not, but contains long-range 
components that cannot be captured with a small size 
cluster Stanescu and Kotliar— . 



proximated in any way using finite range cluster compo- 
nents. 

To chose the proper periodization procedure, we com- 
pare the density of states for cumulant and self-energy 
periodization for both 3 and 4 site clusters. In the vicin- 
ity of half-filling, regardless of the cluster size, the self en- 
ergy periodization results in states lying inside the Mott 
gap, as shown in Fig. [5]anc(6l These states are clearly un- 
physical, as demonstrated by the comparison with the lo- 
cal cluster spectral function, which shows a well-defined, 
clean Mott gap, as shown in Fig. [U On the other hand, 
these mid-gap states are absent in the cumulant peri- 



In contrast, at large dopings in the Fermi liquid phase, 
where both the self-energy and the cumulant are short- 
range quantities, the two methods agree. We conclude 
that the self-energy periodization is appropriate away 
from the Mott transition, while the cumulant scheme 
gives consistent results in a wide range of dopings. In 
the present study we will use the cumulant method re- 
gardless of filling. The reason for the failure of the self- 
energy periodization method is the presence in the half- 
filled regime of self-energy divergences^^ at uo = and 
low temperatures. This divergence of the self energy at 
half-filling is intimately linked to the Mott gap. 
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B. Commensurate Insulators 

Away from half filling, in the electron doped regime, 
there is a substantial discrepancy between the two clus- 
ters. More explicitly, our method predicts that around 
n ~ 1.33 for the triangular and n ~ 1.25 for the rhombic 
cluster, the system becomes an insulator as evidenced by 
the gap in the density of states in Fig. (a) and (b) 
respectively. 

To better understand the commensurate insulator 
phases, we can analyze the resolvents within NCA and 
evaluate the contribution of each cluster state. It is there- 
fore possible to determine which are the dominant chan- 
nels through which the cluster interacts with the bath. 
The relevant quantity is the partial occupancy of a par- 
ticular cluster state 



{flrnm) 



qI 



(0- (14) 



Small occupancy signifies that the corresponding resol- 
vent contributes insignificantly to the cluster spectral 
function and consequently to the hybridization function. 
Therefore, the resolvents with small occupancy can be 
ignored, whereas the ones with large overlap have the 
dominant contribution to the dynamics of the system. 
For the triangular cluster it turns that there are only 7 
resolvents with an overlap larger than 0.01 whereas the 
rest have occupancy less than 0.0025. A similar analysis 
can be performed for the rhombic cluster. As expected 
there are more resolvents with substantial occupancy. In 
Fig. [8] only the dominant ones are plotted which have 
occupancy more than 0.1. Their quantum numbers are 
shown in Table ITTl In all cases, we observe that when the 
dominant resolvents have fillings i/N c close to the lat- 
tice filling, their partial occupancy peaks and transitions 
from and to them become rare, which gives rise to an 
insulating state. In the triangular lattice this has a con- 
sequence when n = 4/3, because in this case there is only 
one dominant resolvent. As a result, there is no apprecia- 
ble overlap with any electronic states, leading thereby to 
a gap in the spectrum. In the rhombic cluster, around the 
critical value n = 5/4 there are two dominant resolvents 
which means that even though the transitions are lim- 
ited, a gap still persists, though it is not as pronounced 
as in the three-site cluster. We have also obtained com- 
mensurate insulating states for fillings n = 2/3 for the 
triangular and n = 3/4 for the rhombic cluster, but in 
this regime convergence is not reached at low tempera- 
ture. The presence of such insulating states ultimately 
points to a limitation of the finite cluster approach to the 
Hubbard model. 

We emphasize that the appearance of the fictitious 
commensurate insulating states is a consequence of us- 
ing small clusters in the numerical calculations. The use 
of these small clusters is dictated by practical reasons. 
However, the interpretation of the results requires spe- 
cial care. Ultimately, the consistency of any result should 
be confirmed by its invariance to cluster size variations. 
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Table I: The quantum numbers of the dominant resolvents for 
the triangular cluster for U=12t, T=0.1t. 
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Table II: The quantum numbers of the dominant resolvents 
for the rhombic cluster for U=12t, T=0.1t. 



Such an example is Mott insulating phase at half filling, 
n = 3/3 = 4/4. The existence and the properties of this 
state can be obtained consistently using various cluster 
sizes. Moreover, as the on-site interaction U is reduced, 
a transition to a metallic state is observed, with some 
small cluster size dependence of the critical parameters. 
In all the present implementations of the CDMFT cluster 
methods for finite dimensions higher than one, there is 
always a difficulty related to the fact that most or all of 
the cluster sites lie on the cluster boundaries. However, 
since large cluster computations are impractical, one has 
to rely on a careful interpretation of the small cluster 
results. We suggest that the ultimate criterion for the 
consistency of small cluster calculations is the indepen- 
dence of the results to the cluster size. 



IV. RESULTS 

In this section we present our main results. In subsec- 
tion IIV At we show that the two-dimensional Hubbard 
model on a triangular lattice is a strongly correlated sys- 
tem by demonstrating that the low-energy physics is con- 
trolled by a weakly dispersing band with spectral weight 
that can be transferred over large energy scales. Sub- 
section IIV Bl is devoted to the characterization of the 
interaction-controlled Mott metal-insulator transition at 
half filling. In contrast to the infinite dimensional case 
when a coherence peak develops inside the Mott gap, we 
find that the insulator-metal transition is characterized 
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Triangle, U=12t, T=0.1t 



Rhomb, U=12t, T=0.1t 
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Figure 8: The occupancies versus rilling for the dominant re- 
solvents in the triangular (a) and rhombic (b) cluster. 
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Figure 9: The evolution of the density of states as a function 
of doping in the upper and lower (inset) Hubbard band. 



by the complete collapse of the Mott gap followed by 
the appearance of a small peak in the density of states. 
Finally, in subsection IIV CI we investigate the system at 
low dopings and discuss a new perspective on pseudogap 
physics. We argue that the pseudogap should not be sim- 
ply identified by the depletion of the density of states at 
the chemical potential, but rather by the change in the 
location of the low-energy modes in momentum space as 
compared with the non-interacting system. According to 
this picture, the pseudogap phase is essentially charac- 
terized by the re-constructed Fermi surface consisting of 
small pockets that vanish in the zero doping limit. 



A. Dispersionless low energy band 

The main figure and the inset in Fig. [9] and Fig. [10] 
reveal a lack of particle-hole symmetry for electron and 
hole doping. This is expected as a triangular lattice does 
not preserve this symmetry. While the asymmetry per- 




2.0 



Figure 10: Evolution of the density of states for the triangular 
cluster around half filling for T — O.lt and U — Ylt in the 
electron doped regime (upper Hubbard band). There is a 
pseudogap feature which does not align with the chemical 
potential (to — 0) and which disappears for a doping of n — 
3.5%. The inset shows the density of states in the hole doped 
regime (lower Hubbard band), which exhibits no pseudogap 
feature. 



sists regardless of the cluster size, the details differ. Of 
particular interest is the presence of a dispersionless sub- 
band residing near the top of the lower Hubbard band 
upon hole doping. The occurrence of such a band is in- 
consistent with Fermi liquid behavior: the chemical po- 
tential crosses the band in an extended area instead of 
at a well defined curve. This band occurs in both the 
triangular and the rhombic cluster, but in the latter it 
appears split. This splitting, shown in Fig. [Til but also 
in the density of states as shown in Fig. [12j may be due 
to the higher resolution gained by the using the rhombic 
cluster. 

A shadow of this band persists unsplit and with less 
spectral weight even in the electron doped regime where 
the system is a normal metal as shown in Fig. [T2j There- 
fore the splitting is due to the fact that it crosses the 
chemical potential. Consequently, in the triangular lat- 
tice, particle hole asymmetry gives rise to a non-Fermi 
liquid behavior for hole doping and metallic behavior in 
the electron doped side. This band structure can be com- 
pared with experimental results Qian et al.— . A full com- 
parison is not possible because only one band, namely the 
a\g is taken into account and the e g > is ignored. However 
the experiment shows the existence of an almost flat band 
with energy — 0.6eV. The present calculation is evidence 
that this band may emerge purely because of strong cor- 
relations. Further evidence that this band arises from 
purely strong electron correlations comes from the fact 
that it is absent in local density approximation (LDA) 
and linear augmented plane wave (LAPW) calculations 
Qian et al. Fig. 3 ofi£. 
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Figure 11: Comparison of the density of states between t 
rhombic and the triangular cluster in the upper and lov 
Hubbard band. There is good agreement in the electr 
doped regime. In the hole doped regime there is a disp 
sionless band which the rhombic clusters resolves as bei 
split. 



B. Mott transition 

The density of states for high enough U exhibits an 
interaction-induced gap around half filling, as shown in 
FigHH Clearly shown in Fig. [14] is the closing of the 
gap between the lower and upper Hubbard bands as the 
on-site interaction decreases, indicative of a Mott transi- 
tion. To pinpoint the precise location of the Mott tran- 
sition, we estimate the Mott gap by the discontinuity, 
A/i, in the chemical potential on either side of half-filling. 
Fig. [15] displays a typical calculation of the chemical po- 
tential as a function of the filling for both 3 and 4-site 
clusters. Because of thermal broadening, this procedure 
would underestimate A/i, which cannot discriminate be- 
low energy scales of the order of kT. For both cluster 
sizes, the results are consistent yielding a gap of A/i w 5t 
for U = 12t. As the inset demonstrates, the discontinu- 
ity in the chemical potential across half-filling vanishes 
at U c ~ 5.7t. However, A/i provides only a rough es- 
timate of the critical U because the precise magnitude 
of the gap is obscured as shown in Fig. [15l To probe 
the transition more directly, we plot the density of states 
for different values of U around the estimate obtained 
from the chemical potential analysis. Fig. [16] displays 
clearly that for U > U c , with U c w 5.75t and U c w 4.75£ 
for the triangle and rhombic cluster respectively, the sys- 
tem is an insulator, whereas for U < U c a Drude peak 
emerges at the chemical potential. Note that the density 
of states near the chemical potential remains unchanged 
(relative to its value for U > U c ) although the density 
of states at the chemical potential develops a non-zero 
value as a transition is made to a metallic state. Conse- 
quently, the states that fill in the Mott gap and give rise 
to the coherence peak arise from spectral weight transfer 
from high energy, the essence of Mottness. This is illus- 
trated in Fig. [13 which compares the spectral function 
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Figure 12: The spectral function as a function of k and to for 
U = 12£, T = O.lt and n = 0.99. Upper panel: triangular 
cluster, lower panel: rhombic cluster. Both methods give a 
dispersionless band, but in the rhombic cluster it is resolved 
in two. 



right before and after the transition. Such a redistribu- 
tion of spectral weight far from the chemical potential 
near the Mott transition has been recently observed in 
the manganites^. 

The operative mechanism for the Mott transition in 
the three and four-site cluster analysis stands in contrast 
to the scenario predicted by DMFT 12 . In this scenario a 
coherent peak of constant height exists at the chemical 
potential, which successively narrows as U increases to 
U c . For U > U c , the peak vanishes and the upper and 
lower Hubbard bands become well separated. CDMFT 
in the plaquette offers a different scenario^: first a pseu- 
dogap opens at the chemical potential which smoothly 
grows to form a full Mott gap at the critical U. In the 
triangle, there is no formation of a pseudogap. Instead 
the coherence peak at the chemical potential loses weight 
Z as U increases and is smoothly replaced by a gap which 
broadens as the two bands separate. 



C. The pseudogap 

In light of the physics in the cuprates, one of the main 
questions that needs to be addressed is whether or not 
a Mott system on a triangular lattice exhibits a pseu- 
dogap. We focus here on the single-particle density of 
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Figure 15: The chemical potential as a function of filling for 
the 3 and 4 site cluster at U = 12t and T = O.lt. There 
is agreement anywhere apart from the commensurate fillings 
n — 4/3 and n — 5/4 respectively. Inset: the evolution of 
the Mott gap as a function of filling for the triangular cluster. 
The gap is linear in U and vanishes at U c = 5.74t. 



Figure 13: Spectral function for U = 12t, T = O.lt and 
n — 1.04 for the triangular (upper panel) and rhombic cluster 
(lower panel). The chemical potentials crosses a narrow band 
and the behavior is metallic. The dispersionless sub-band at 
the top of the lower Hubbard band persists, but with less 
spectral weight and unsplit.. 
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Figure 14: The density of states in the vicinity of half filling 
for Triangular cluster. The gap closes with U. 
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Figure 16: (a) density of states around the Mott transition for 
the triangular cluster at n — 1 T = O.lt. As U decreases, the 
Drude peak around the chemical potential gradually looses 
weight until it turns turns into a gap. The transition occurs 
at some U c ~ 5.75t. (b): density of states around the Mott 
transition for the rhombic cluster at T = 0.12t and n — 1 and 
U c w 4.5t. 



states as a function of filling. In the triangular lattice, 
there is only indirect evidence from a boson analysis of 
the optical conductivity for a pseudogap in the paramag- 
netic phase (n < 1.5)^^ which disappears in the strange 
metal phase (n > 1.5). 

Our results for the single-particle density of states on 
either side of half-filling are summarized in Figs. [9] and 
fTUl Although a dip-like features exists for both the 3 and 
4-site clusters, it is displaced from the chemical potential. 
For higher dopings, the density of states is smooth in the 
vicinity of the chemical potential. Consequently, we find 
an absence of a pseudogap near half-filling on a triangular 
lattice. This result is consistent with the Quantum Monte 
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Figure 17: The spectral function at half filling and T = O.lt 
as a function of to and k along the r — > K — > M — > T path for 
U = 6t (upper panel) and £7 = 5.5t (lower panel) as obtained 
from the triangular cluster. The spectral weight mostly from 
energy around 4.5t is transferred to the chemical during the 
transition. 
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Figure 18: Evolution of the DOS around the chemical poten- 
tial for n — 1.02 using the triangular cluster. The pseudogap 
feature is dispersed as the temperature increases. 



Carlo calculations Kyung 10 as the impurity solver cou- 
pled with CDMFT on the triangular lattice. We traced 
the origin of the dip-like feature to an energy splitting 
of two resolvents in the triangular lattice (three in the 
rhombic) with total spin differing by 1. The relevant en- 
ergy scale is J = 6t 2 /U. To probe this feature further, 
we calculated its evolution as the temperature is lowered. 
As is evident from Fig. [18j as the temperature is lowered, 
the dip-like feature remains below the chemical potential 
and more importantly the peak at the chemical poten- 
tial sharpens, as would be expected for a metallic state. 
Consequently, using the dip-like feature in the density of 
states as the criterion for the existence of a pseudogap, 
we conclude that no pseudogap exists for the Hubbard 
model on a triangular lattice for either electron or hole 
doping. 

However, this analysis is incomplete. Let us approach 
the pseudogap problem from a different perspective. In 
the undoped cuprates, the quasiparticle dispersion be- 
low the Mott gap is characterized by four maxima at the 
(±7r/2, ±7r/2) points in the Brillouin zone, as revealed 
by ARPES measurements on Ca 2 Cu02Cl^. This fea- 
ture is also present in the half-filled Hubbard model on 
a square lattice^. For a weakly (hole) doped system, 
if one adopts the naive picture of a rigid band shift, 
one would expect the chemical potential to move to the 
top of the lower Hubbard band and intersect it some- 
where in the vicinity of the four minima. The resulting 
Fermi surface would consist of four small hole pockets in 
the vicinity of the (±7r/2, ±7r/2) points, while the rest 
of the large Fermi surface observed in optimally doped 
and overdoped cuprates (or at large doping values in 
the calculations for the Hubbard model) would be com- 
pletely obliterated. That is, the low-energy excitations 
are gaped everywhere in the Brillouin zone, except on 
the boundary of the small Fermi pockets. This picture 
seems to be consistent with ARPES measurements on 
underdoped cuprates^, as well as the infrared Hall ef- 
fect Rigal et al.— , Shi et al.— and quantum oscillation 
measurements^. Nonetheless, strong coupling calcula- 
tions show that the naive rigid band picture is, in fact, 
incorrect and that strong correlations play a crucial role 
in pseudogap physics. 

One of the essential aspects of strong correlations is 
spectral weight transfer. To illustrate its role, let us ap- 
proach the pseudogap problem for Hubbard model on 
a triangular lattice starting from the Mott insulating 
phase. Note that, in contrast to the square lattice model, 
in this case particle-hole symmetry is always absent and 
the antiferromagnetic interactions are frustrated. In Fig. 
([T9|) we show the top of the lower Hubbard band (left 
panel) and the bottom of the upper Hubbard band (right 
panel) for a half-filled system with U = I2t. The two 
bands are separated by a Mott gap of about 5£, and 
the chemical potential sits in the middle of the gap. 
The two surfaces are defined by the smallest frequencies 
at which the spectral function exceeds a certain small 
threshold 8 A. Variations of 5 A produce only small shifts 
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of the two surfaces, but their shapes remain essentially 
the same. First, let us focus on the upper Hubbard band. 
As shown in Fig. ([T9|h it is characterized by a set of min- 
ima along a large closed curve around the T point, not 
far from the Fermi surface corresponding to that of the 
non-interacting half-filled system. A small electron dop- 
ing would move the chemical potential near this set of 
minima of the upper Hubbard band. In the rigid band 
picture, one ends up with two almost circular Fermi sur- 
faces that define a narrow electron ring. However, our 
strong-coupling calculation leads to a different picture. 
Shown in Fig. (|2T|) is the spectral function along the 
r — > K — > M — > T path in the Brillouin zone and a small 
frequency window about the chemical potential for two 
small doping values. For comparison, we also show the 
bottom of the upper Hubbard band for the insulator (up- 
per panel) within a similar frequency window. Note that 
for the insulator the bottom of the band corresponds to 
some middle points between T and K and between M 
and T. In the vicinity of the T point, there is no spec- 
tral weight within our frequency window [see also Fig. 
(|T9|) ]. However, once we dope the system, some spectral 
weight is transferred from high energies to low energies, 
so that a low energy band clearly forms in the vicin- 
ity of T just below the chemical potential. As a result, 
the strongly correlated narrow band that controls the 
low-energy physics disperses across the chemical poten- 
tial generating a large Fermi surface consistent with the 
Luttinger theorem. A crucial difference from the square- 
lattice case is that the anisotropy along this large Fermi 
surface (or along the line defining the minima of the up- 
per Hubbard band for the insulator) is very weak. By 
contrast, for the square lattice there is a qualitative dif- 
ference between the (7r/2,7r/2) and the (0, n) regions of 
the Brillouin zone as illustrated, for example, by the ex- 
istence of minima near the nodal points. This lack of 
anisotropy leads to the sudden appearance of a large 
Fermi surface upon doping and thus to the absence of 
a pseudogap. 

The situation appears somehow different in the case 
of the lower Hubbard band. As shown in Fig. ([T9|h the 
top of the lower Hubbard band is extremely flat and ex- 
tends over a significant portion of momentum space, all 
around the boundary of the Brillouin zone. Very weak 
maxima can be identified near the K points. In this case, 
the rigid band picture would suggest that a weakly hole- 
doped system is characterized by low-energy excitations 
that extend over a large portion of the Brillouin zone and 
that small Fermi pockets would possibly form near the 
K points and extend rapidly with doping. Again, in the 
vicinity of T there is no spectral weight at low energy. 
However, in contrast to the upper Hubbard band, this 
lack of low-energy excitations at small momenta persists 
upon doping. Shown in Fig. [20] is the low-energy spectral 
function along the same r — > K — > M — > T path for the 
insulator (upper panel) and two values of doping. The 
relevant spectral weight transfer contributes this time to 
the re-shaping of the low-energy narrow band that exists 



at momenta far from the V points. At increased doping 
values, this band becomes more dispersive and generates 
a large Fermi surface consistent with the non-interacting 
Fermi surface of a system with the same filling factor. 
Nonetheless, at very small doping values, several ques- 
tions remain. First, it seems that the chemical potential 
crosses the narrow band in an extended area of the Bril- 
louin zone, rather than along a well-defined Fermi line. 
This stands in sharp contradiction with Fermi liquid the- 
ory. However, one has to take into account that our re- 
sults are obtained at a finite temperature of order O.lt, 
and thus the energy resolution is severely limited. To 
establish exactly the position of the Fermi surface at low 
dopings would require a much better energy resolution, 
and consequently a much lower temperature, would be 
necessary. The second question concerns the existence of 
a pseudogap. One typically understands the pseudogap 
as a reduction in the number of low-energy modes below 
a certain energy scale. Here we propose a slightly differ- 
ent view. We define the pseudogap phase as a physical 
state occurring close to a Mott insulating state and char- 
acterized by the existence of small Fermi pockets with an 
area proportional to the doping level x = 1 — n. By con- 
trast, a normal Fermi liquid is characterized by a large 
Fermi surface with an area that is related, via the Lut- 
tinger theorem, to the filling n. Consequently, one should 
view a system in the pseudogap phase as a doped Mott 
insulator. At the same time, the system represents a re- 
normalized Fermi liquid characterized by a reconstructed 
Fermi surface. From this perspective, the weakly hole- 
doped Hubbard model on a triangular lattice is in the 
pseudogap state. Within the energy and momentum res- 
olution of the present method, the Fermi surface appears 
as a set of small pockets around the K points that expand 
rapidly upon doping. A normal Fermi liquid is estab- 
lished at a doping level of a few percent. We emphasize 
that a crucial condition for the realization of this pseu- 
dogap phase was the existence of the small anisotropy in 
the lowest energy excitations of the Mott insulator. To 
study in detail the formation and the evolution of the hole 
pockets, calculations using larger clusters (i.e., having a 
better momentum resolution) and lower temperatures are 
necessary. 



V. CONCLUSIONS 

We have studied a strongly correlated electron sys- 
tem on a triangular lattice using an implementation of 
the NCA+CDMFT scheme. The key technical aspects 
of this implementation are presented in detail. Numer- 
ical results are obtained for two types of clusters con- 
taining three and four sites, respectively. We stress that 
the cluster size analysis is a required step in any cluster 
DMFT-type calculation and argue that the relative in- 
variance of the result against increasing the cluster size 
is the ultimate consistency criterion. The fundamental is- 
sue concerns the short-versus long-range character of the 
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n=1.00 






Figure 20: Comparison of the spectral function at half fill- 
ing (top), 1% (middle) and 3%(bottom) hole doping, for the 
rhombic cluster at U = 12£, T = O.lt. 




Figure 21: Comparison of the spectral function for half filling 
(top), 1% (middle) and 3%(bottom) electron doping, for the 
rhombic cluster at U = 12t, T = O.lt. 



electron correlations and the nature of the quantity that 
properly describes them. We find that the self-energy is 
not a short-range quantity in the vicinity of half-filling 
and therefore cannot be captured using the cluster com- 
ponents. However, within our momentum and energy 
resolution, we find that the cumulant satisfies the local- 
ity requirements and can be used for re-constructing the 
lattice quantities. In this context, a very high priority 
for future cluster DMFT studies should be to establish 
the relevant range for the self-energy and the cumulant 
in various parameter regimes. Larger cluster calculations 
are required to clarify this point. Nonetheless, the task 
is of pivotal importance because, if in a certain regime, 
both the self-energy and the cumulant are long-ranged 
quantities, the presently available real- and momentum- 
space cluster DMFT schemes are not applicable. 
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At low doping values, we find that the Hubbard model 
on a triangular lattice is strongly correlated with low- 
energy physics controlled by a quasi-dispersionless band. 
As a result of correlations, the band is very narrow, and 
its spectral weight can be transferred over large energy 
scales. A band with such features cannot be described 
by any weakly coupled approach. We also find that a 
metal-insulator transition occurs at a critical value of the 
on-site interaction U c ~ 5.6 ± 0.15t, which depends very 
weakly on the size of the cluster. This value is much lower 
than the critical interaction determined of U c ~ 10. bt in 
CDMFT calculations using exact diagonalization as the 
impurity solver Kyung 10 , but it is closer to U c « 6.9t, 
which is the critical interaction obtained by continuous- 
time Monte Carlo Lee et al.— . Finally, we discussed 
the pseudogap problem in the context of the Hubbard 
model on a triangular lattice. In contrast to the square- 
lattice case, we find no evidence for a dip in the density of 
states positioned at the chemical potential. However, a 
momentum-resolved analysis shows that the locus of the 
low-energy excitations of the weakly hole-doped system 
is qualitatively different from that of a non-interacting 
system. Therefore, we propose a new framework for 
discussing the pseudogap phenomenon, which in essence 
involves a momentum-dependent characterization of the 
low-energy physics, rather than a momentum-integrated 



one. We define a pseudogap state as a state characterized 
by low-energy excitations occurring only in a relatively 
small region in momentum space, qualitatively different 
from the location of the low-energy quasiparticles of the 
non-interacting system, and having an area that shrinks 
to zero when approaching the Mott insulator. Conse- 
quently, the system in the pseudogap state is charac- 
terized by a re-constructed Fermi surface consisting of 
small pockets. We find that the conditions necessary for 
the appearance of these pockets is a strongly momentum- 
dependent self-energy which produces quasiparticles with 
anisotropic properties along the Fermi surface. There- 
fore, the pseudogap is intrinsically linked to Mott physics 
as emphasized recently , which is the source of the long- 
range self-energy. Within the resolution of the present 
calculation, we find that the momentum-dependence of 
the self-energy is much weaker for the triangular lattice, 
as compared to the square lattice, leading to a pseudogap 
only in the very weak hole-doped regime. 
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